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ABSTRACT 

We generalize the Marinari-Parisi definition for pure two dimensional quantum gravity (k = 
2) to all non unitary minimal multicritical points (k > 3). The resulting interacting Fermi gas 



! theory is treated in the collective field framework. Making use of the fact that the matrices 

^ ' evolve in Langevin time, the Jacobian from matrix coordinates to collective modes is similar to 

^ ■ the corresponding Jacobian in d = 1 matrix models. The collective field theory is analyzed in 
OS ! 

. the planar limit. The saddle point eigenvalue distribution is the one that defines the original 

i \ multicritical point and therefore exhibits the appropriate scaling behaviour. Some comments on 

J!h \ the nonperturbative properties of the collective field theory as well as comments on the Virasoro 

■ constraints associated with the loop equations are made at the end of this letter. There we 
3 ■ 

%j . also make some remarks on the fermionic formulation of the model and its integrability, as a 
nonlocal version of the non linear Schrodinger model. 
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Hermitian one matrix models defined at a single pointQ describe at their double scaling 
multicritical points minimal nonunitary conformal matter coupled to two dimensional quantum 
gravity § @ §. 

Pure gravity (k = 2) as well as all higher multicritical points of even order k (k = 4, 6, ...) 
exhibit non perturbative ambiguities, instabilities and violations of the Schwinger- Dyson (loop) 
equations § § § QQ. 

These problems can be traced back to the fact that the critical matrix potentials of even 
order k are bottomless. Thus a sensible definition of these models should correspond to well- 
defined stabilized matrix potentials. 

Bottomless matrix potentials occur also in certain multimatrix models describing unitary 
matter coupled to two dimensional quantum gravity such as the two matrix model corresponding 
to the Ising case [[UJ . Thus, the problem of stabilization is associated not only with one matrix 
models. 



Marinari and Parisi |TT| have suggested a possible way out of this difficulty in the case of 
pure gravity by supersymmetrizing the model. 



The bosonic sector hamiltonian in [|4]] is also the forward Fokker-Planck hamiltonian asso- 
ciated with the Langevin equation whose force term is minus the gradient of the original zero 
dimensional matrix action. 



Therefore, the definition of pure two dimensional quantum gravity in [ITJ is equivalent to 
the stabilization procedure developed in [15 for bottomless actions, as far as the bosonic sector 
of the former is concerned. 



The stabilized pure gravity model of [[LI]] was further analyzed by |13|] , where the one eigen- 
value double scaled hamiltonian was extracted, and a nonperturbative ambiguity free analysis 
of the density of particles and density of state was made. 

The forward Fokker-Planck hamiltonian^], used in JTTJ] [[0|] for the pure gravity case, reads 

3 i.e. - matrix models defined over a zero dimensional space time. We will refer to them as zero dimensional 
models. The stabilized models, in which matrices depend upon the Langevin time coordinate will be referred 
to as the one dimensional models. One should not, of course, mix this terminology with "dimension" counting 
in target space - i.e. values of the central charge of matter coupled to gravity. 

4 In [Q the existence of such an TL was postulated, without specifying its details to define "time or- 
dered" correlators in the zero dimensional theory 
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for a general matrix potential U(<&) 

H = \ T ^-m + l u '^w + f u ' m] (1) 

Here $ is an hermitian matrix, depending on the Langevin time coordinate (i.e. the bosonic 
coordinate of superspace used in |TT|). 

7i in eq. (1) is a well defined Schrodinger operator. Its potential is clearly bounded from 
below and grows to plus infinity as the matrix eigenvalues become infinite. 

Therefore, a well defined unique normalizable ground state vector ^ r o(^ > ) exists. Indeed, if 
is bounded from below such that the Boltzman weight of the zero dimensional matrix 
model is normalizable this ground state is given by 

tf ($) = J_ e ^T r Um ? z= f d JV* $e -0ZW(*) (2) 
y Z J 

In this case the vacuum energy is strictly zero and supersymmetry is not broken. Moreover, 
expectation values of operators, all at infinite Langevin time project only onto the ground state 
^o, and are identical to the corresponding correlators in the original zero dimensional matrix 
model: 

< # | X ($) - • ■ O n ($) | * >= \ J d Jva $e-WOi($) ■ ■ ■ n ($) (3) 

If U ($) is unbounded from below, the zero dimensional Boltzman weight is unnormalizable 
and the corresponding matrix model exists only at a saddle point level. Supersymmetry is 
broken and the vacuum energy Eq is positive. Alternatively - the appropriate Langevin equation 
has only runaway solutions and the Fokker-Planck probability density at any finite portion of 
matrix eigenvalue space, decays asymptotically in Langevin time t as e~ Eot ]0| |ft6HR. 

However, Fokker-Planck averages of operators normalized by the Fokker-Planck average of 
the unit operator are well defined as the Langevin time goes to infinity, and correspond to 

<0l($) •••0n($)>*~oo= /A|$ ($)| J dW-Wu (4) 

5 For this property to hold also in the double scaling limit, we must ensure that the vacuum remains nonde- 
generate even as N — > oo, i.e. - that the energy eigenvalue E\ of the first excited state of Tt does not coalesce 
with Eq. As was shown in |fl| the mass gap E\ — Eq double scales in the WKB approximation for k = 2. 
Moreover, we will show that it double scales for any value of k. Thus the vacuum state remains nondegenerate. 
This was tacitly assumed in |l2| ]. 
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Here ^o^) is the normalizable ground state of 7i, and all the operators on the r.h.s. of Eq. 
(4) are at t — oo. We thus consider Eq. (4) as the stabilized definition for correlators in case 
of bottomless matrix potentials. 

As is well known, the laplacian over hermitean matrices acquires the form 

_ Tr _^_ = L_ y A(x) + ( nonsinglet U (N) angular \ 

<9$ 2 A(x) ~^ dxf \ momentum terms / 

where Xj are the matrix eigenvalues and A(xi) is the Vandermonde determinant. This leads 
to the mapping of eigenvalue dynamics onto that of a one dimensional Fermi gas [ 18| . Clearly, 
the ground state \Po(<I>) mentioned above is a U(N) singlet. 

For a generic potential U(<&), the hamiltonian in Eq. (1) contains a two body interaction 
term among eigenvalues^ 

4 <9$ 4 Xi - Xj 

Thus, generally - the one dimensional gas of eigenvalues is an interacting Fermi gas. 

The matrix potential of minimal degree that leads to the k = 2 multicritical (pure gravity) 
point is U2(<&) = — 1$ 3 + |$ 2 + $ where we have followed the normalizations of |19| . The critical 
coupling constant corresponding to the k = 2 point is A c = 1. The matrix potential t^t^) is 
bottomless. Therefore the stabilized theory will exhibit spontaneous supersymmetry breaking. 



This issue was analyzed in pl |. In this case, the interaction term in Eq. (6) reduces to an 
interaction of the eigenvalues with a constant background field, proportional to the number of 
eigenvalues, N, namely, Hint = 4 T r (A$ — 1). 

Therefore, the singlet sector of Eq. (1) reduces effectively to a system of N non interacting 
Fermions in one dimensionF] [TT| [Tj| 

N 1 02 N 

n singlet = 1 E(-^2 7T2 + V{xi)) = (3 2 Y,h(x l ,p i ) (7) 



6 Note that for £/'($) = J2o< n C n Tr$ n we obtain H mt = £ <„ C»E£o Tr^Tr^™" 1 ^, therefore from 

the point of view of the (noncritical) one dimensional matrix theory Ho = ^Tr(— -J^ + ^-t/'($) 2 ), whose 
eigenvalues form a noninteracting Fermi gas in the singlet sector, TLi n t may be interpreted as higher curvature 
terms |20| that push the system to its multicritical point. 

7 In N/ (3 was set to 1 and A was varied, while in Q A was set to its critical value A c = 1 and N/ (3 was 
varied around its critical value (N/(3) c = 1. 
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where the external one body potential is 

V{x)= l -{U' 2 {x)f + \ j(Xx-l) (8) 

The ground state of the Fermi gas described in Eqs. (7) and (8) is obtained by filling the first 
iV one particle levels of V(x). The associated Fermi energy must be evaluated self-consistently 
from the N-dependent potential V(x). Therefore, unlike the case of the d=l model, the Fermi 
energy is not a free parameter that can be used to define the double scaling limit. 

In order to study the higher stabilized multicritical points (k > 3), one has to cope with the 
interaction term in Eq. (6). Since only the U(N)-singlet ground state ^o^) of 7i is involved - 
it is natural to analyze the interacting gas in terms of the Fermion density operator - i.e., the 
collective field <p(x) associated with the matrix <3>. 

Following p2[ |23| we define f\ the collective field as 

1 1 N 

4>{x) = - T r 5(l . x -§) = -Y J 5{x-x i ) (9) 

P P i=i 

which implies that <f>(x) is a non negative operator and obeys the normalization condition 

N 



r iv 

J <t>{x)dx = - (10) 



P 

Since $ depends on the Langevin time t and its dynamics is fixed by Ji in Eq. (1), the 
Jacobian of the transformation from the matrix eigenvalue variables to the collective field is 
exactly the same Jacobian as in the d = 1 case |23] f\ 

In this letter, we concentrate on the planar approximation to Ji in order to establish the 
fact that the k — th order multicritical behaviour is respected by our formalism. 

Discussing the exact loop equations and the associated non perturbative effects is deferred 
to a subsequent publication 



8 We have used the normalization of |^6| . 

9 The use of precisely this Jacobian is dictated by the very definition of the supersymmetrized model in 
accordance with [^4| |25|] , bypassing the need to invoke arguments of the type used by (24j] , or postulates about 
the form of the zero dimensional partition function as in j25| , that are needed if one makes the transformation 
to collective modes directly in the zero dimensional matrix model. 
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The planar collective field action So[<f>] = /3 3 J C^dxdt for the matrix Hamiltonian TC in 
Eq. (1) is given by [|2| |23 . 



y ,s, x ,U ' (x) N(/ N1 lr, U'(x) - U'iy) , , 



:n) 



Here fip is a Lagrange multiplier (the chemical potential) that enforces the constraint of Eq. 
(10). Unlike the case of d = 1 matter [p^j , ^ is not a free parameter whose deviation from a 
critical value is used to define the scaling behaviour. Eq. (11) yields the planar contribution 
to the genus expansion. The term ^(U'(x)) 2 is the external effective potential in which the 
eigenvalues move. Its contribution to the action clearly produces the ^-T r ({7'($)) 2 term in Eq. 
(1). Similarly, the non local bilinear interaction term in Eq. (11) is the collective field version 
for Tiint in Eq. (6). The cubic term in Eq. (11) arises due to the transformation from matrix 
eigenvalues to collective modes. If one considers the first few terms in the 1/(3 expansion of 
the collective field action, including the kinetic term that is identical to the one used in d = 1 
matrix models, one obtains (up to ambiguities known in this expansion) a non local collective 



field theory analogous to the one developed in [23 



However, unlike the latter, fluctuations of our collective field theory around the WKB 
solution do not correspond to genuine string field components due to the fact that the target 
space dimension is negative. Nevertheless, it might include minor fractions of string field 
compontents - namely, discrete states. 

The planar collective field equation of motion is readily found to be 

6S[<f>] tt 2 n2 U\xf 



S6(x) i 

v 7 \ planar 



Pl y dy u'(x)-u'(y) 



<j>(y) = (12) 



x-y 

where 4>(x) is subjected to the constraint of Eq. (10), and that by definition, 4>(x) is nonnegative. 

A crucial observation is that for % = 1 and fip = 0, this nonlinear nontrivial integral 
equation is identical to the planar limit of the Schwinger- Dyson equation obeyed by the loop 
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operator in the original zero dimensional matrix model [I9| Jl8| Jl7j 
F(z) 2 - U\z)F{z) + rj(z) = 

F(z)= lim <±T r ^—>, V (z) = lim < -T r U ' {z) - ^ ($) > 



(13) 



when z approaches the real axis. 

The fact that Eq. (12) is identical to the planar loop equation of the original matrix model 



is not surprizing and conforms with the postulates of JT2[. Moreover, it seems that the WKB 
expansion of Eqs. (3) and (4) should correspond term by term to the genus expansion of the 
corresponding Schwinger-Dyson (i.e. loop) equations in the original modef^] 

Therefore, under the conditions ^ = 1 and fip = 0, (b(x) that solves Eq. (12) is just the 
planar limit eigenvalue density of the original Dyson gas in an external potential U(x). 

Thus, for a matrix potential U(Q) in the universality class of the k-th multicritical point, 
4>(x) will exhibit k-th order multicritical behaviour. In particular, if we chose [[HJQ 



U' k {x) = ^y^[(2 - *) fc - V - m + (14) 

then <f)(x) is supported only along the closed segment [-2,2] on the real axis and is given by 

^ ) = (t^(4-^(2-x)- J^2 (15) 



and satisfies 



2 

4>k{x)dx = 1 . (16) 

-2 
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As z = x - ie , e -> 0+ we have fh9 [h 8l 1' 



F(z) = ^U'(x) +iiv<f)(x) 

thus making the imaginary part of Eq. (13) vanishing and its real part proportional to Eq. (12). 

n It can be shown [^Tj, for example, that an Ehrenfcst theorem, associated with Eq. (4), given by 
< alhi ^(z^*) ah >= leads order by order in the WKB expansion (in the N < (3 phase) to the genus 
expansion of the corresponding loop equation in the original matrix model, of which Eq. (13) is the leading 
(planar) term. 

12 The (+) subscript means that in an expansion of the r.h.s. of Eq. (14) around x — oo we keep only 
nonnegative powers of x. 
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Eqs. (14)-(16) are the solution to Eqs. (10) and (12) precisely at the k-th multicritical 
point. 

In order to have a complete solution of the stabilized model on the sphere one has to 
show that these equations are properly deformed by turning on the cosmo logical constant. 
Namely, that under a deformation of TrUk{§) by the pucture operator Tr$, Tr[4($) — > 
Trf/fc($) + fisTr®, where /i# is the bare cosmological constant, there exists a solution to Eqs. 
(10) and (12), supported along a single segment on the real eigenvalue axis, such that 1 — ^ 
scales as p- 2k / 2k + 1 . Such a deformation of the normalization condition in Eq. (16) may be 
obtained by allowing the multicritical end-point (i.e. x = 2) of the support of (f>k(x) in Eq. (15) 
to vary, on a proper scale.0 This scale must be that of the double scaled fluctuations of the 
matrix near its critical point $ c = 2, i.e., [@] (3~ 2 / 2k+l . 

Therefore, the required eigenvalue distribution should be supported along a segment [-2,b] 
where 

6 = 2- e(3~ 2/2k+1 (17) 

in which e is a finite real parameter. 

Such a deformation of <fik{x) alone is not enough to obtain the desired scaling behaviour of 
1 — since it generically induces all k-1 relevant deformations [[I]] JTP| present at the k th multi- 



critical point. The desired solution to Eqs. (10) and (12) must therefore include counterterms 



that will cancel these unwanted scaling contributions to 1 — Thus, it must have the general 



form 



0( x ) = 9l{2 + x )h{b-x) k ~^ + Y^(3 Z ^ 1 af»^(2 + ^(6-x)^ 



^ — j TT 



on its support [-2,b]. 

Here C n = ^jfeffi normalizes f^ 2 4>n(x) to unity when b — > 2 as in Eqs. (15) and (16). 

2(k-n) 

j3 2k + 1 (n < k — 1) are the scaling dimensions of the relevant perturbations at the k-th 
multicritical point |1|] [^] ensuring that <fi(x) scales as a whole as j3~ 2k + 1 when b—x ~ y(3^ 2 ^ 2k+1 . 

13 In the d=l matrix model ^3|, variation of the chemical potential /if changes the location of the classical 
turning points which are the end points of swpp{(j)(x)} . 



8 



Finally, af£> (n < k — 1) are the double scaled couplings of the relevant scaling operators, 
that will be uniquely fixed by the scaling behaviour of 1 — 4. 

Indeed, using the elementary integral J^ 2 ^(2 + x)^{b — x) n ~^dx = (^) n+1 Eqs. (10), 
(17) and (18) yield 

j = (i - ^ k+1 + e «i fc) A fc -(i - -\r +i (19) 



2 



where we have set A = (3 2k + 1 . 



In the vicinity of the k-th multicritical point, 4- scales as 

N 2k 

— = l - tf3~—i (20) 

where t is the renormalized cosmological constant. Therefore, expanding the r.h.s. of Eq. (19) 
in powers of A, the coefficients of A, A 2 , • • ■ A fc-1 must vanish, providing k — 1 equations that 
uniquely fix the k — 1 unknowns (a n ) in Eq. (18). The latter are readily found to be given in 
closed form by 

Jk) _ f k + 1 \ d\k-n 



whence 4 = 1 — [k + l)(|) fc j3~ 2k+1 ■ (1 + 0{f3^ 2k + 1 )). Therefore, the renormalized cosmoligical 
constant is given by0 t — (k + l)(|) fc , or equivalently 

e = 4(-^ T ) 1/fc (22) 

Eqs. (17), (18), (21) and (22) comprise together the desired deformed solution 4>(x) to Eq. 
(12). 

To this eigenvalue distribution we may add relevant scaling deformations with arbitrary 
couplings, that will not alter the normalization condition in Eq. (10) fL9| . 

Eq. (22) exhibits the well known scaling relations QXJ] on the sphere between the specific heat 

e of the original matrix model and the cosmoligical constant t. That is, Eq. (22) implies that 

the string susceptibility is given by 7 str = —1/k. Equivalently, it implies that it is impossible 
14 Eqs. (21) and (22) were obtained in closed form by M. Moshe. 
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to obtain a negative cosmological constant in our solution for even values of k, keeping e real 
- i.e., that it holds only in the N < (3 phase of the theory. 

General solutions of the planar loop equation (Eq. (12)) in the N > (3 phase of multicritical 
points of even order k would probably have branch point singularities in the complex eigenvalue 
plane (arranged in complex conjugate pairs) other than the two branch points on the real axis 
H ||, giving rise to multicut (planar) macroscopic loops which may be supported at negative 
loop lengths where they are also oscillatory and nonpositive definite). The latter fact, should 
however, in our opinion, be considered as a sickness of the planar loop equation, or equivalently 
of the saddle point (WKB) solution of the stabilized matrix model - rather than as an a priori 
signal of non-perturbative instabilities of the Marinari-Parisi procedure. 

f3(f)(x) given by Eqs. (18) and (21) clearly equals the WKB approximation to the particle 
density of the interacting Fermi gas of eigenvalues (integrated up to the Fermi energy). 

2k — 1 2 

By construction, <f>(x) scales as /3~ 2fc + 1 in terms of the scaling variable y = (b — x). 
Therefore, according to well known arguments, the WKB density of one-particle states at the 
Fermi energy, corresponding to one-particle excitations of the ground state of the gas is given 
by0 

l,dN. r b dx ^2fe- 3 r°° dy n ik-a , . . . 

I ~ (32E+T / y = f3—^p(e) (23) 



p\ dE )\K=B W J- 2 <p(x,b) Jo <f>(y,e) 

This integral is generically well behaved, and diverges only as e — >• 0, i.e., as one approaches 
the /c-th multicritical point .0 

Thus, the energy-gap of these excitations double-scales as 



Af e ) = [-(—), l" 1 = /ri+T — ~| 1-— \^ f'M) 



2fc-3 



where the last proportionality is expected on general scaling arguments. For odd values of k, 
our solution <f)(x) is well defined either in the N < (3 phase or in the N > (3 phase. It has the 
same singularity structure either for positive e and t or for negative ones. Thus, in such cases 

A(e) in Eq. (24) is well defined in both phases. 

15 Such one-particle states should exist at least in a frame work of Hartree-Fock analysis of the gas. 
16 (j)(x) given by Eqs. (18) and (21) obviously vanishes only at x = b, —2, where it has generically Wigner 
semicircle singularities as long as e 7^ 0. 
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In case of even values of k, <j>{x) we have found exists only in the N < (3 phase where A(e) 
is real and positive. Recall that a positive real double scaled excitation gap is essential for the 
stabilization of the model by Eqs. (3) and (4) to work. 

For k = 2 our results reproduce the WKB analysis of JTTj] |TJj || in the N < (3 phase of the 
matrix model. Indeed, in this case, Eqs. (17), (18), (21) , (22) and (24) read 

6(x) = { i [(2 + x) (2 - e/?" 2/5 -xp(2 + §/T 2/5 -x) -2 < x < 2 - e /T 2/5 m) 
VK ' \ otherwise 1 ; 

as well as f = 1 - 3(f) 2 /?" 475 and A(e) = $ p-W = ^(| f )^-i/5 „ (1 _ |)|. 

This single segment supported <j>(x) corresponds to the fact that the Fermi energy of the 
eigenvalue gas described by Eq. (7), coincides with the shallow local minimum of V(x) in Eq. 
(8) near x = +2. This eigenvalue distribution leads to a single cut maxroscopic (planar) loop 
operator |1| whose Laplace transform is supported only at non-negative loop lengths |§ |28| . 



Analysis of Eqs. (7) and (8) in the N > (3 phase was made in |J. We have briefly described 
it in our comments following Eq. (22), where we have also stated our interpretation of the 
difficulties pointed out in [[| fj\ 

We have remarked above 7 that the analysis of the k = 2 point in |H| were performed by 
keeping the original cubic matrix potential critical while varying N/ f3 around its critical value 
1. At this point of our discussion it is clear that []13| has been successful in doing so because 
precisely for that cubic potential does N//3 couple to its appropriate scaling perturbation in 
Hi n t ~ namely, the puncture operator Tr$. 

Up to this point we have established that all stabilized multicritical one matrix models 
are equivalent on the sphere to the corresponding original (zero dimensional) models. This is 
expected also to hold to all orders in the genus expansion |TJ| |^7| (at least in the N < [3 phase 
in the case of even k). This conclusion is a good starting point and a motivation to study the 
non-perturbative nature of the stabilized multicritical models. 

It is well known that the set of all multicritical points of (zero dimensional) one matrix 



models form an integrable hierarchy - namely that of the KdV equation |1] |2| |29|] [fL4}| 



17 We have also constructed two solutions to Eq. (12) for k — 2 and ^ > 1 in which pip ^ 0. One of these 
solutions with a single segment support exhibits scaling behaviour of the k = 3 point. The other solution with 
two real and two complex conjugate branch points seems to correspond to the one discussed in . 
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The latter may be represented by a set of Virasoro constraints that annihilate the exact 
partition function of the matrix model |j30| . 

Thus, in order to establish quantitative results concerning nonperturbative differences be- 
tween the original and stabilized matrix models one must first try to reformulate the latter as 
an integrable hierarchy 

Such a construction, if possible, may be carried either by attempting to formulate the loop 
equations of the stabilized models in a manner analogous to or by checking whether the 
interacting Fermi gas itself forms an integrable system. 

In the first case, one makes an explicit use of the basic observables of the matrix model - 
namely the loop operators - hence an immediate comparison of the two types of models may 
be done. It may well be that using the exact collective field theory in this respect, without 



expanding its Jacobian in powers of 1/(3, turns out to be quite valuable |25[, especially due 
to the fact that this Jacobian (for the stabilized model) is identical to the one used in d=l 
matrix models |24j|. The enormous symmetry possessed by the collective field theory in the 
latter case |31| or by its Fermionic counterpart |32| might have counterparts in our case as well 
(even though surely in our case we have more complicated Fermi "Fluid dynamics"). 

The other method proposed above will yield, if it turns out to be successful, the entire 
spectrum of the interacting Fermi gas and its exact S-matrix. This may be by itself an interest- 
ing result in the theory of 1 + 1 dimensional integrable models - extending the local nonlinear 
Schrodinger model into a nonlocal version (with special one and two body interactions - derived 
from multicritical potentials)^ 



Summary and Conclusions 



We have shown that the Marinari-Parisi definition of pure gravity (k=2) may be extended to 

stabilize all higher one matrix multicritical points. This was done by demonstrating the equiva- 

18 Note also that by quantizing the interacting Fermi gas - in a Hartree-Fock formalism one obtains, expanding 
around the WKB (i.e., planar) Fermi field configuration, a non local version of the the Thirring model, in an 
analogous manner to the Dirac theory obtained in [p6{ . 
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lence on the sphere of the original and stabilized models. It seems that this equivalence should 
hold to all orders in the genus expansion. The collective field of the stabilized matrix model 
turned out to be useful in coping with eigenvalue interactions in our semiclassical treatment of 
the Fermi gas. 

We have pointed out that, since the whole structure of multicritical points of one matrix 
models may be transferred to the Marinari-Parisi arena, the most important questions are 
whether they form there an integrable hierarchy as the original models do, and if the latter is 
answered on the affirmative - how does it differ from the original KdV hierarchy? 
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